Stable, Strongly Attractive, Two-State Mixture of Lithium Fermions in an Optical Trap 
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We use an all-optical trap to confine a strongly attractive two-state mixture of lithium fermions. By 
measuring the rate of evaporation from the trap, we determine the effective elastic scattering cross 
section Ana 2 to show that the magnitude of the scattering length \a\ is very large, in agreement with 
predictions. We show that the mixture is stable against inelastic decay provided that a small bias 
magnetic field is applied. For this system, the s-wave interaction is widely tunable at low magnetic 
field, and can be turned on and off rapidly via a Raman tv pulse. Hence, this mixture is well suited 
for fundamental studies of an interacting Fermi gas. 
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Trapped, ultracold atomic vapors offer exciting new 
opportunities for fundamental studies of an interacting 
Fermi gas in which the temperature, density and interac- 
tion strength can be independently controlled. Recently, 
a degenerate gas of fermionic 40 K has been produced by 
using a two-state mixture to enable s-wave scattering and 
evaporation in a magnetic trap jl|. By removing one 
species, the properties of the noninteracting degenerate 
gas were measured, demonstrating that the momentum 
distribution and the total energy obey Fermi-Dirac statis- 
tics |l[. However, the properties of interacting two-state 
fermionic vapors have not been explored experimentally. 

Theoretical treatments of an interacting Fermi gas 
have focused extensively on 6 Li [g-||]. Certain two- 
state 6 Li mixtures are predicted to be strongly attractive, 
i.e., they have anomalously large and negative scattering 
lengths pi arising from a near-zero energy resonance in 
the triplet state pl| . It has been predicted that these 
strongly attractive mixtures can undergo a transition to 
a superfluid state at a relatively high transition tempera- 
ture [ppl]. In addition, the two-state effective interaction 
potential is widely tunable in a magnetic field, permit- 
ting systematic studies of fundamental phenomena such 
as collective oscillations for both the normal and super- 
fluid phases Pj5|,|| , as well as new tests of superconduc- 
tivity theory W. 

Unfortunately, magnetically trappable mixtures in 6 Li 
with large s-wave scattering lengths are not stable, since 
there are correspondingly large spin-exchange and dipo- 
lar decay rates pMpJ. Hence, the methods employed to 
study degenerate K are not applicable. For this reason, 
we developed an ultrastable CO2 laser trap to confine a 
stable mixture of the two lowest 6 Li hyperfine states p2| . 
However, attaining a large and negative scattering length 
in this mixture requires high magnetic fields B > 800 G 
to exploit either a Feshbach resonance or the triplet scat- 
tering length @,0. 

In this Letter, we show that there exists another sta- 
ble hyperfine state mixture in 6 Li which has the following 
unique properties. First, we predict that the scattering 
length a is large, negative, and widely tunable at low 
magnetic field B. By monitoring the rate of evaporation 



from the CO2 laser trap at a fixed well depth, we measure 



\a\ = 540±ll Q a at B = 8.3 G. This result confirms for 
the first time that very large scattering lengths exist in 
6 Li mixtures. The predicted scattering length is — 490 ao 
at B — 8.3 G, consistent with our observations, and is 
expected to increase to — 1615 ao as B — » 0. Second, 
we find that this system is stable against spin exchange 
collisions provided that B^0. In addition, the dipolar 
decay rate is predicted to be very small Q, consistent 
with our observations. Finally, in the experiments, a Ra- 
man 7r pulse is employed to abruptly create an interacting 
mixture from a noninteracting one, a desirable feature for 
studies of many-body quantum dynamics I14J . 
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FIG. 1. 6 Li hyperfine states, labeled |1) to |6) in order of 
increasing energy in a magnetic field. The magnetic quantum 
number of each state is denoted by m. The hyperfine constant 
a hf = 152.1 MHz. 

Fig. p] shows the hyperfine states for 6 Li labelled |1) — 
1 6), in order of increasing energy in a magnetic field. At 
low field, the states |1) and |2) correspond to the \F — 
1/2, m) states, while states |3) through |6) correspond to 
states \F = 3/2, m). At nonzero magnetic field, only the 
magnetic quantum number m is conserved. The subject 
of this paper is the |3) — |1) mixture. 
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FIG. 2. Magnetic field dependence of the scattering length 
ti3i for a mixture of the |3) and |1) hyperfine states of 6 Li. 

Fig. g shows the scattering length 031 for the |3) — |1) 
mixture as a function of magnetic bias field B. We es- 
timate (131(B) by using the asymptotic boundary con- 
dition (ABC) approximation fin] . This calculation in- 
corporates the singlet and triplet scattering lengths [Oj , 
as = 45.5±2.5 ao and ax = — 2160±250 qq, and a bound- 
ary radius which we take to be R = 40 ao |10| . The scat- 
tering length varies from — 1620 ao (~ 3ot/4 as B — > 0) 
to — 480 a at B = 10 G. The results of our approximate 
calculation for B = to B = 200 G are confirmed within 
10% by van Abeelen and Verhaar jl|] using a coupled 
channel calculation which includes the uncertainties in 
the potentials. At higher fields, near 800 G, we believe 
the scattering length exhibits a Feshbach resonance (not 
shown). Above this resonance, the scattering length ap- 
proaches the triplet scattering length of — 2160 ao. 

The 1 3) — |1) mixture is stable against spin-exchange 
collisions provided that a small bias magnetic field is 
applied. Spin-exchange inelastic collisions conserve the 
two-particle total magnetic quantum number My, where 
Mt = —1 for the |{3, 1}) state. Note that {,} denotes 
the antisymmetric two-particle spin state, as required 
for s-wave scattering which dominates at low tempera- 
tures. There are no lower-lying antisymmetric states with 
Mt = — 1. Hence, exothermic collisions are precluded. 
The only other states with Mt = — 1 are |{4, 2}) and 
I {5, 3}). Without an adequate bias magnetic field, transi- 
tions to these states lead to population in level |4). Then, 
exothermic |{3,4}) -> |{3,2}) and |{4,1}> -> |{1,2}) 
collisions can occur. With an adequate bias magnetic 
field, the energy of states |{4, 2}) and |{5,3}) can be in- 
creased relative to that of state |{3,1}) by more than 
the maximum relative kinetic energy, i.e., twice the well 
depth during evaporative cooling. By energy conserva- 
tion, spin-exchange transfer is then suppressed. In this 
case, the inelastic rate is limited to magnetic dipole- 
dipole (dipolar) interactions which contain a rank 2 rela- 
tive coordinate operator of even parity J2J . Since parity is 
conserved and p — wave — > p — wave scattering is frozen 



out at low temperature, the dominant dipolar process is 
a small s -> d rate in which |{3, 1}) -> |{1, 2}) |f]. 

In the experiments, the CO2 laser trap is initially 
loaded from a magneto-optical trap (MOT) fl2[| . At 
the end of the loading period, the MOT laser beams are 
tuned near resonance and the intensity is lowered to de- 
crease the temperature. Then, optical pumping is used 
to empty the F = 3/2 state to produce a 50-50 mixture 
of the |1) — |2) states. These states are noninteracting 
at low magnetic field, i.e, the scattering amplitude van- 
ishes as a result of an accidental cancellation (To| . With 
a CO2 laser trap depth of 330 //K, up to 4 x 10 5 atoms are 
confined in the lowest-lying hyperfine states at an initial 
temperature between 100 and 200 fiK. A bias magnetic 
field of 8.3 G is applied to split the two-particle energy 
states by ~ 16 MHz. This is twice the maximum attain- 
able energy at our largest well depth of 400 ^fK = 8 MHz. 
After a delay of 0.5 second relative to the loading phase, 
a pair of optical fields is pulsed on to induce a Raman 
it pulse. This pulse transfers the population in state |2) 
to state 1 3) in two microseconds, initiating evaporative 
cooling in the resulting |3) — |1) mixture. The optical 
fields are detuned from resonance with the D2 transition 
by ~ 700 MHz to suppress optical pumping. If the Ra- 
man pulse is not applied, the trapped atoms remain in 
the noninteracting |1) — |2) mixture and exhibit a purely 
exponential decay with a time constant ~ 300 seconds. 

An acousto-optic modulator (A/O) in front of the CO2 
laser controls the laser intensity, which is reduced to yield 
a shallow trap depth of 100 /jK. By using a shallow well, 
we avoid the problem that the elastic cross section be- 
comes independent of the scattering length at high en- 
ergy, as described below. In addition, the shallow well 
greatly reduces the number of loaded atoms and makes 
the sample optically thin, simplifying calibration of the 
number of trapped atoms. To determine the trap pa- 
rameters, the laser power is modulated and parametric 
resonances [ |15[ are observed at drive frequencies of 2v for 
three different trap oscillation frequencies v. At 100 [iK 
well depth, v x = 2.4 kHz v y = 1.8 kHz and v z = 100 Hz, 
where the trap laser beam propagates along z. Using the 
measured total power as a constraint, we obtain the trap 
intensity 1/e 2 radii, w x = 50 /im and w y — 67 /im, and 
the axial intensity 1/e 2 length, Zf ~ 1.13 mm, where Zf 
is consistent with the expected Rayleigh length within 
15%. 

The number of atoms in the trap N(t) is estimated 
using a calibrated photomultiplier. The detection sys- 
tem monitors the fluorescence induced by pulsed, retrore- 
flected, <7± probe and repumper beams which are strongly 
saturating (I/I sa t = 26 for the strongest transition). 
To simplify calibration, only the isotropic component of 
the fluorescence angular distribution is measured: The 
collecting lens is placed at the magic angle |16| of 55° 
(P2 (cos 0) = 0) with respect to the propagation direction 
of the probe beams. The net efficiency of the detection 
system is determined using laser light of known power. 
The primary uncertainty in the calibration arises from 



the excited state population fraction, which we estimate 
lies between 1/4 and 1/2. 
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FIG. 3. Number of trapped atoms versus time for evapo- 
ration of a 1 3) — |1) mixture of 6 Li at a fixed well depth of 
100 fiK. The solid curve shows the s-wave Boltzmann equa- 
tion fit for a scattering length of |asi| = 540 ao- Inset: 0-1 
second. 

Fig. [3] shows the number of trapped atoms N(t) mea- 
sured for the |3) — 11) mixture at a well depth Uq = 100 [iK 
and a bias field of 8.3 G as a function of time between 
5 ms and 20 seconds after evaporation is initiated. For 
times beyond 50 seconds (not shown), the evaporation 
stagnates, and we observe an exponential decay of the 
cooled 1 3) — |1) mixture with a time constant of 370 sec- 
onds over a period of a few hundred seconds. The error 
bars are the standard deviation of the mean of ten com- 
plete runs through the entire time sequence. 

A model based on the s-wave Boltzmann equation p"7J ] 
is used to predict N(t) for comparison to the experiments. 
This equation is modified to include the density of states 
for a gaussian potential well [12[ and to include the en- 
ergy dependence of the elastic cross section. Assuming 
a short range potential and a symmetric (s-wave) spatial 
state, the cross section takes the form 



a(k) = 



87ra§ 1 



„2„2 • 

„ a 31 



(1) 



where Tik is the relative momentum. For fc|a3i| << 1, 
the cross section is maximized. When fc|a3i| >> 1, the 
cross section approaches the unitarity limit 8ir/k 2 which 
is independent of 031. Note that fc|a3i| = 1 corresponds 
to a relative kinetic energy of e = Ti /(2/i a^), where /i = 
M/2 is the reduced mass. For |o,3i | = 500 ao, e = 115 /iK. 
For a two-state mixture of fermions, the effective cross 
section is reduced from that of Eq. [j] by a factor of 2 since 
pairs of colliding atoms are in an antisymmetric hyper- 
fine state with a probability 1/2. This effective cross 
section is used in a Boltzmann collision integral for each 
state i = 1, 3. A decay term —Ni(t)/r with r = 370 sec 



is added to account for the measured trap lifetime. A 
detailed description of our coupled Boltzmann equation 
model will be published elsewhere. 

The coupled s-wave Boltzmann equations for the two 
states are numerically integrated to determine N(t) using 
the well parameters as fixed inputs. From the calibrated 
photomultiplicr signal, assuming that 1/3 of the atoms 
are in the excited state, we obtain an initial total number 
Nq = 44, 000. For this case, the initial collision rate 
in Hz is estimated to be 1/(27tt c ) ~ NqM oqv 3 / '(fc^T) , 
where v 3 — v x v y v z , <tq = 8770,3^, and M is the 6 Li mass. 
Assuming la.31 1 = 500 a , r c = 30 ms. Hence, for t > 0.3 
seconds, when on average 10 collisions have occurred, the 
sample should be thermalized as assumed in the theory. 

The best fit to the data starting with 22,000 atoms in 
each state is shown as the solid curve in Fig. ft The \ 2 
per degree of freedom for this fit is 1.4 and is found to be 
very sensitive to the initial temperature To of the atoms 
in the optical trap. ^From the fit, we find To = 46/zK, 
which is less than the well depth. We believe that this 
low temperature is a consequence of the MOT gradient 
magnet, which is turned off after the MOT laser beams. 
The effective well depth of the optical trap is therefore 
reduced until the gradient is fully off, allowing hotter 
atoms to escape before the Raman pulse is applied to 
create the |3) — |1) mixture. The fit is most sensitive to 
data for i > 0.5 second, where the thermal approxima- 
tion is expected to be valid. From the fit, we obtain the 
scattering length la.31 1 = 540±25ao, which is within 10% 
of the predictions of Fig. Q The quoted error corresponds 
to a change of 1 in the total x 2 - 

We determine the systematic errors in 031 due to the 
uncertainties in the calibration and in the population im- 
balance as follows. The data is fit for an initial num- 
ber of atoms Nq of 58,000 and 29,000, corresponding to 
an excited state fraction of 1/4 and 1/2. This yields 
|a3i I = 440 ± 20 ao and la.31 1 = 750 ± 42 ao, respectively. 
Note that for the larger scattering lengths, the cross sec- 
tion given by Eq. |lj approaches the unitarity limit and 
the error increases. We assume that the initial popula- 
tion imbalance for states |3) and |1) is comparable to that 
of states |2) and |1) in the optically pumped MOT. To 
estimate the latter population imbalance, we use state- 
selective Raman 7r pulses to excite |2) — > |3) or |1) — v |6) 
transitions in the MOT. Probe-induced fluorescence sig- 
nals from states |3) or |6) show that the initial |1) and 
1 2) populations are equal within 10%. Note that resid- 
ual population in state |2) is expected to be stable and 
weakly interacting, since we estimate (032! < 30 ao for 
< B < 50 G using the ABC method, and a 12 ~ [0. 
Using the parameters for the fit shown in Fig. @, but 
changing the initial mixture from 50-50 to 60-40, we find 
a slight increase in the fitted scattering length from 540 ao 
to 563 ao. Thus, the uncertainty in the calibration of the 
number of atoms produces the dominant uncertainty and 
|a 3 i| = 540± 2 oo , a . 

To demonstrate that evaporative cooling is occuring, 
rather than just trap loss, we have also measured the 



final temperature of the mixture using release and recap- 
ture [Hi from the C0 2 laser trap. We obtain 9.8 ± 1 fiK, 
which is within 10 % of the final temperature of 8.7 /xK 
predicted by the Boltzmann equation model. An excel- 
lent fit to the data is obtained for the final temperature, 
which describes a thermal distribution. However, the ini- 
tial temperature is not so readily measured, as it is non- 
thermal before evaporation is initiated, and is rapidly 
changing during evaporation, unlike the final tempera- 
ture, which stagnates. 

Good fits to the evaporation data are obtained neglect- 
ing inelastic collisions, suggesting that the dipolar rate 
for the 1 3) — |1) mixture is small, in contrast to the scat- 
tering length. A limit on the dipolar loss rate for the 
|3) — |1) mixture can be estimated from the r = 370 sec- 
ond lifetime of the mixture after evaporation stagnates. 
For equal populations in both states, dipolar decay re- 
sults in an initial loss rate h = — Gn 2 /4, where G is the 
dipolar rate constant and n is the total density. To ob- 
tain a high density, the trap is loaded at a well depth 
of 330 /iK and the temperature of the atoms is reduced 
by evaporation to T ~ 30 ± 1 /iK. The number of atoms 
remaining in each state after evaporation is estimated to 
be N = 6.5±2.2 x 10 4 , where the uncertainty is in the 
calibration. We cannot rule out the possibility that one 
state is depleted on a long time scale, since we do not 
directly measure the individual state populations. How- 
ever, we believe that, after evaporation stagnates in the 
deep well, a |3) — |1) mixture remains, since subsequent 
reduction of the well depth yields final temperatures con- 
sistent with evaporative cooling. Note that the mixture 
ratio is not critical: An 80-20 mixture yields an initial 
loss rate h = —0.16 Gn 2 , ~ 2/3 that of a 50/50 mixture. 
For a fixed 330 /iK trap depth, v 3 = 2.6 ± 0.3 kHz 3 , and 
the phase space density for one state in the harmonic 
approximation is then p p h = N/(kBT/hv) 3 — 7 x 10~ 4 . 
This corresponds to a maximum total density of n = 
2p ph /\ B = 6.4 x 10 n /cm 3 , where \ B = h/y/2irM k B T. 
Since the exponential decay time of the |3) — |1) mixture 
is similar to that obtained in the noninteracting |1) — |2) 
mixture, we assume the loss is dominated by background 
gas collisions. Thus, we must have Gn/A « 1/r, which 
yields G « 2 x 10 _14 cm 3 /sec. This result is consis- 
tent with the value G ~ 2 x 10 -15 cm 3 /sec predicted for 
the dipolar rate constant at 30 /iK by van Abeelen and 
Verhaar |l3| . 

Future experiments will employ continuous evapora- 
tion by slowly reducing the well depth jLf|. In this case, 
very large scattering lengths can be obtained at low tem- 
peratures and small well depths by using a reduced bias 
magnetic field B. By adiabatically recompressing the 
well, experiments can be carried out with the precooled 
atoms in a deep trap to obtain high density as well. In 
such experiments, the final low temperature limits the 
number of atoms in the high energy tail of the energy dis- 
tribution, exponentially suppressing spin-exchange colli- 
sions for B 7^ 0. For example, if a total of 3 x 10 5 atoms 
were contained in our trap at a well depth of 400 /iK, the 



Fermi temperature Tp — 7/iK and the Fermi density is 
4 x 10 13 /cm 3 . At a temperature of T = 0.1 T F = 0.7 /iK, 
a bias field of B = 0.16 G would split the two-particle 
hyperfine states by ksTp + 12 fc^T, suppressing the spin 
exchange rate by exp(— 12), and giving 031 ~ — 1200 ao 
Alternatively, as shown in Fig. g, large 031 can be ob- 
tained at moderate B ~ 300 G. 

In conclusion, we have observed that an optically 
trapped |3) — |1) mixture of 6 Li atoms has a very large 
scattering length at low magnetic field. This mixture 
is stable against spin-exchange collisions provided that 
a small bias magnetic field is applied. The evaporation 
curves measured for this mixture are in good agreement 
with a model based on an s-wave Boltzmann equation 
which neglects inelastic processes. We have predicted 
that the scattering interactions are strongly attractive 
and widely tunable at low magnetic field. If the param- 
eters described above for deep wells can be attained, the 
system will be close to the threshold for superfluidity g] 
and ideal for investigating frequency shifts and damp- 
ing in collective oscillations J3|Jq]. Further, since s-wave 
interactions can be turned on and off in a few microsec- 
onds, this system is well suited for studies of many-body 
quantum dynamics. 
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